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Abstract 

Using data collected by LHCb experiment in 2011 we report the measurements 
of charge asymmetries in charmless decays of B mesons in two or three charged 
kaons or pions. We find positive charge asymmetries in the channels B® — > 
K~tt + (3.3a), B + ->■ K+TT+iT- (2.8a), and B + ->■ vr+vr+vr- (4.2a) and negative 
in B° -> K+vr- {6a),B + -»• K + K + K~ (3.7a), B + -»• K + K~ir + (3.0a). 

1 Introduction 

We present recent studies from the LHCb experiment on direct CP violation in B 
meson charmless decays. Using data collected in 2011 we have studied the decays: 
B° -»■ K+7T- and B° s -> K~n + [I]; B + ->■ K+tt+tt- and 5+ -»■ K+K + K~ g] 
and i? + K + K~ir + and 5 + — > 7r + 7r + 7r~ [3] (complex-conjugate modes are implied 
except in asymmetry definitions). The main observable in all the analyses is the charge 

asymmetry in decay rates into a final state /, defined as A C p = f?§~fefrf~^ • These 
channels have been extensively studied by the B factories and the Tevatron [3] but 
only now, with the large samples collected by LHCb, it was possible to establish first 
evidence (significance > 3a) of CP asymmetry in four of the channels. The exceptions 
are the well established B° — > K + n~ [1] for which we report on the most precise 
measurement and the B + K + tt + 7t~ for which our measurement for the inclusive 
Acp has a significance of 2.8a. The main features of the LHCb experiment that made 
this possible are: the very efficient particle identification and vertex location system, 
dedicated hadronic trigger and the ability to reverse the magnet polarity on the data 
taken [S]. 

CP violating transitions and possible new physics effects originating in loop dia- 
grams are accessible in all channels. The decays B + — > K + tt + tt~ , B + — > K + K + K~ 
and B° — > K + n~ proceed through similar CP conserving penguin diagrams and have 
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access to CKM CP violating transitions b — > u at tree level. In these cases the 
smallness of penguin amplitudes is compensated by ~ A 2 factor as compared to the 
factor ~ A 4 for the tree contribution, where A ~ 0.22 is the expansion factor of the 
Wolfenstein parametrization of the CKM matrix. CP violation effects can potentially 
be large. The decays B + — > K + K~tt + , B + — > 7r + 7r + 7r~ and B° s — > K~ir + access CP 
violating transitions b — > u and b — > t in tree and penguin diagrams respectively, both 
at the same order ~ A 3 leading to the overall penguin contribution much smaller than 
the tree and observation of CP violation is disfavored. 

Using 0.35 fb _1 of pp collisions at = 7 TeV, LHCb established the first evidence 
of CP violation in the B® system in its decays to Kir pairs, and also reported the most 
precise measurement of the CP asymmetry of B° — > Kir. These results are presented 
in Section 2. 

The three-body decays offer a good environment to study CP violation enriched 
by the interference pattern of the two-body resonances in the Dalitz plot. The charge 
asymmetry in regions of the phase space is likely to be greater than the integrated 
ones. With the complete data set of 1 fb _1 collected in 2011, we measure inclusive 
CP asymmetries for the three-body channels as well as their distributions in the 
phase space. It is a consequence of CPT invariance and unitarity that the sum 
of all partial decay widths of a particle should match that of its antiparticle. In 
fact, CPT is more restrictive. If we divide all possible decays of a particle into 
classes containing final states that are mutually distinct under the strong interactions, 
CPT imposes the particle-antiparticle match should be observed in each class [6]. 
Motivated by the physics, by CPT invariance that connects the decays K ± n + n^ to 
K ± K + K~ and K + K~tt ± to n^^^ final states in two distinct classes, and also 
due to the experimental resemblance, it was natural to group the three-body analysis 
in two sets: i? ± — > K ± 7t + ti^ and 5 ± — > K ± K + K~ are described in Section 3 and 
B ± ->■ K + K~tt ± and B ± ->■ tt ± tt + 'k~ in Section 4. 

2 Direct CP violation in B° s -> K~7T + and B° -> 
K + tt~ decays 

The probability for a b quark to decay as B® — > Kir is about 14 times smaller than 
that to B° — > Kir, and a stronger rejection of the combinatorial background was 
required. The polarity of the LHCb magnetic field is reversed from time to time 
in order to partially cancel the effects of instrumental charge asymmetries. Of the 
0.35 fb _1 , about 43% was acquired with one polarity and the rest with opposite 
polarity. The LHCb hadronic trigger is particularly efficient for the two-body modes. 
It requires one very high transverse momentum good quality track in the event. To 
distinguish the Ktt final states we rely on the RICH particle identification system. The 
efficiencies and misidentification rates were carefully controlled using large calibration 



2 



data samples of D* — > Dix —t (Ktt)ott and Ab — > pir decays. The raw asymmetry is 
computed from maximum likelihood fits to the Ktt invariant mass spectra shown in 
Figure [TJ The raw asymmetry needs to be corrected for two residual effects. The first 
is induced by the detector acceptance, event reconstruction and the strong interaction 
of the final state particles with the detector material. The second is associated with 
the production of the B mesons and takes into account the dilution factor due to 
mixing. The overall correction to the raw asymmetry is AAcp{B°) = — 0.007 ± 0.007 
and AAcp(B^) = 0.010 ± 0.002 , where the errors are statistical. The final measured 
CP asymmetries are 

A CP (B° S Kir) = 0.27 ± 0.08(stat) ± 0.02(syst) , 
A C p(B° Kir) = -0.088 ± 0.011 (stat) ± 0.008(syst) , 

(1) 

where the systematic uncertainties are dominated by the correction effects and in the 
case of the B® — > Kir there is also a leading contribution from the combinatorial 
background description. 




Figure 1: Invariant Ktt mass spectra obtained using the event selection adopted for best 
sensitivity to (a),(b)A C p (B° ->• Ktt) and (c),(d)A C p (-B° Ktt). Plots (a) and (c) 
represent the K + tt~ invariant mass whereas plots (b) and (d) represent the K~ tt + invariant 
mass. The results of unbinned maximum likelihood fits are overlaid. The main components 
contributing to the fit model are also shown. 
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3 Direct CP violation in B + ->> K + tt + 7t~ and B + — >> 
K + K + K~ decays 

With relatively large branching fractions of (9(10~ 5 ) the B + — > K + tt + tt~ and B + — > 
K + K + K~ signals have comparable data samples and the analysis differs only in 
the particle identification and in specific peaking background estimation. The odd 
number of kaons in the final state and kinematics similarity enable a straight forward 
use of the large B ± — > J/tpK ± control sample to account for residual effects from 
the production mechanism, as well as the differences in the K + and K~ interaction 
with the detector material. The raw charge asymmetry, A^p W , calculated from the 
positive and negative event yields and the corrected A C p is given by A C p(K ± h + h~) = 
A^ w (K ± h + h-)-A^ w (J/i/jK ± ) + A CP (J/^K ± ) where hh stands for KK or mr and 
the current PDG value for Acp{J/ijj K^) is (1 ± 7) x 10~ 3 [I]. The subtraction was 
done separately in sub samples of the data divided by trigger selections, in order 
to cancel a possible trigger asymmetry for kaons. The total event yields are listed 
in Table 1 determined from unbinned extended maximum likelihood fits to the B^ 
candidate invariant masses in the range 5079 — 5479 MeV/c 2 . 

The final results for the inclusive charge asymmetries of B^ 1 — > K ± 7r + 7T~ and 
B ± ^K+K- are: 

A CP (B ± K ± ti + ^-) = +0.034 ± 0.009(stat) ± 0.004(syst) ± 0.007(7/^^) , 
A CP (B ± ^K+K-) = -0.046 ± 0.009(stat) ± 0.005(syst) ± 0.007( J/V^) , 

where the first uncertainty is statistical, the second is systematic, and the third is the 
uncertainty of the A CP of B ± J/ipK^ from the PDG gj. The main contributions to 
the systematic uncertainty are attributed to the use of the control channel to account 
for residual asymmetries to signal model and to the acceptance in the Dalitz plot. 

The two relevant variables used to represent the Dalitz plot are the two-body in- 
variant masses (m^ +7r _, m|- ±7rT ) in the case of B ± — )■ K ± 7r + n~ decay and {jn 2 K + K - low , 
m 2 K + K - high ) in the case of B ± — > K ± K + K~ decay, where the symmetrical phase space 
is folded. To visualise the asymmetries in phase space, the Dalitz plots of candidates 
with three-body invariant masses within ±40 MeV/ c 2 of the peak were divided into 
bins with equal population using an adaptive binning algorithm. In each bin, the 
variable Aq P = was computed from the number of B ± event candidates, N ± , 



mode 




K ± K + K~ 




7r ± 7r + 7r 




B- 
B+ 


18 168 ± 170 
17540 ± 169 


10 289 ± 110 

11 606 ± 117 


30 140 ± 179 
30 984 ± 182 


2718 ± 71 
2111 ±66 


619 ±47 
875 ± 50 



Table 1: Event yields for B and B + samples and for the whole data set. 
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Figure 2: .A^p in Dalitz plot bins with equal event populations for — > K ± K + K (top 
left) and B ± K ± 7r + 7r (top right). The asymmetry Aq P is calculated from the bin con- 
tents including backgrounds, and is not corrected for production or detection asymmetries. 
^raw £ or m 2 K+K _ j (bottom left) and m^ +?r _ (bottom right) projections. 



which includes signal and background events in that bin and is not corrected for effi- 
ciency. The resulting A^ P distributions over the Dalitz plots are shown in the upper 
plots of Figure |2j Large asymmetries are clearly visible in the low m 2 K+K ^ low and low 
m^ +7r _ regions respectively for B ± — >■ K ± K + K~ and B ± K ± -n + it~ . 

To analyse the CP asymmetries over the phase space without background influence, 
the raw asymmetries were plotted in bins of the two-body invariant mass projections, 
shown in the bottom plots of Figure [2j We observe large opposite sign asymmetries 
in the low m 2 r+7T - and m 2 K+K _ low invariant masses not clearly associated to a resonant 
state. 



4 Direct CP violation in B + K + K~tv + and B + ->> 
7r + 7r + 7r~ decays 

The samples of B + — > K + K~n + and B + — > tt + tt + it~ are of comparable sizes, much 
smaller than the B + — > A' + 7r + 7r~ and B + — > K + K + K~ ones and with much higher 
background levels. We measure the raw asymmetry defined in terms of the B + and 
B~ event yields and corrected by the acceptance in the Dalitz plot as ACC A^p W = 
N M±, where R is the ratio between the B + and B~ data- weighted average 



(N B _/R)+N B+ 

efficiency in the Dalitz plot, found to be R nnn = 1.023 ± 0.029 and Rkk-k = 0.991 ± 
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0.027. The total signal yields are listed in Table 1. 

Here the net strangeness in the final state is zero and the residual detection asym- 
metry refers to a pion and is measured by LHCb to be = (0.00 ± 0.25)% [7] 
using large D*^ calibration samples. The residual production asymmetry is esti- 
mated from the B ± — > J/tpK ± control channel, eliminating the factor relative to 
the detection associated to bachelor kaon measured from calibration samples, to be 

= -0.010 ±0.002 [TJ. 

In order to cancel out a possible trigger asymmetry for hadrons, the final asymme- 
tries are calculated from data divided into two subsamples depending on the trigger 
selection. The inclusive CP asymmetries are the weighted averages of the trigger 
subsamples and are measured to be 

A CP (B ± -+ TrW") = +0.120 ± 0.020(stat) ± 0.019(syst) ± O.OO7(J/0.K' ± ) , 
A CP (B ± -> K+K^) = -0.153 ± 0.046(stat) ± 0.019(syst) ± 0.007( J/xpK*) , 

where the first uncertainty is statistical, the second is the experimental systematic 
uncertainty, and the third is the uncertainty of the Acp value of — > J/ipIC* 1 from 
the PDG 0]. 

Investigating the distribution of asymmetries in the Dalitz plots we observe large 
negative contributions at low m 2 K+K _ invariant mass as well as large positive asym- 
metries at low m 2 +7r _, particularly in the high fn 2 r + JT - high region, as can be seen in 
Figures |3] and |H To quantify this effect we compute the local CP asymmetry in 
the phase space regions in the same way as the inclusive one, correcting for the ac- 
ceptance, production asymmetry and detection asymmetry. The CP asymmetry of 
B ± — > 7T ± 7r + 7r _ in the region m 2 +7r _ low < 0.4GeV 2 /c 4 and "i 2 ^- high > 15GeV 2 /c 4 is 
determined to be 

A CP (B ± -+ 7rV"7r" region) = +0.622 ± 0.075(stat) ± 0.032(syst) ± 0. 007 (J/ipK*) . 

The CP asymmetry of B^ — >• K + K~tt ± in the region m 2 K+K _ < 1.5GeV 2 /c 4 is mea- 
sured to be 

A CP (B ± -> K+R-n* region) = -0.671 ± 0.067(stat) ± 0.028(syst) ± 0.007(7/^^) . 

5 Conclusion 

We have observed a rich pattern of CP violation in the charmless decays of the 
B mesons. Contrary to naive expectations, we observe large CP violation in the 
B + — > K + K~7i + , B + — > 7r + 7r + 7r~. Intriguingly large opposite sign Aqp is observed 
in the low K + K~ and invariant mass regions of the three-body decays. The 

effect seems not to be associated to any resonat state. LHCb will triple its sample 
with the data collected in 2012. 
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Figure 3: I? 1 * 1 — >■ K + K^ir ± event yields (signal and background) as a function of 
m^- +x _ (left), where the empty triangles represent B~ and the filled triangles repre- 
sent B + ; fits of B~ (center) and B + (right) for candidates with m 2 K+K ^ < 1.5 GeV 2 /c 4 . 
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Figure 4: B ± — )• 
for m 2 + 



7r ± 7r + 7r 
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event yields (signal and background) as a function of 
high > 15GeV 2 /c 4 (left), where the empty triangles represent B~ 



and the filled triangles represent B + ; and fits of B (center) and B 
candidates with m 2 + _ low < 0.4GeV 2 /c 4 and m 2 r+7r _ u h > 15 GeV 2 /c 4 . 
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